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' We have mapped two regions at the northern tip of the Magellanic Stream in neutral 

Ci ■ hydrogen 21-cm emission using the Arecibo telescope. The new data are used to study 

J> . the morphology and properties of the Stream far away from the Magellanic Clouds, as 

^ , well as to provide indirect constraints on the properties of the Galactic Halo. We inves- 

H ' tigate confinement mechanisms for the Stream clouds and conclude that these clouds 

■ ■ ■ ' cannot be gravitationally confined or in free expansion. The most likely mechanism 

for cloud confinement is pressure support from the hot Galactic Halo gas. This al- 
lows us to place an upper limit on the Halo density: nh(15 kpc) = 10"'^ cm^^ and/or 
rzh(45 kpc) = 3 X 10~^ cm~^ depending on the distance. These values are significantly 
higher than predicted for an isothermal stratified Halo. 
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ABSTRACT 



1. Introduction 

The Magellanic Stream (MS) is a thin, ^ 10° wide (Putman & Gibson 1999), tail of neutral 
hydrogen (HI), emanating from the Magellanic Clouds (MCs: the Large Magellanic Cloud, LMC, 
and the Small Magellanic Cloud, SMC) and trailing away for almost 100° across the sky (— 60°< 
Dec < +15°). This huge HI structure is the most fascinating signature of the wild past interaction 
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of our Galaxy with the MCs, and the MCs with each other. At the same time, streaming so far 
away from the MCs into the Galactic Halo, the MS is an excellent probe of the properties of the 
outer Galactic Halo, for which otherwise we have very few observational clues. 

For many years the MS was viewed as a complex of six discrete concentrations (Mathewson 
et al. 1974). The new HI Parkes All-Sky Survey (Putman & Gibson 1999) reveals the more complex 
nature of the MS, with a fascinating network of filaments and clumps. No stars have been found so 
far in the MS, making it the closest medium to a primordial environment that exists in the Local 
Group today. 

Ha emission appears to be routinely detected throughout the MS (Weiner et al. 2002), suggest- 
ing that the MS is being ionized by photons escaping from the Galaxy (Bland-Hawthorn & Maloney 
1999). While this mechanism is in accord with the Ha observations of many high velocity clouds 
(HVCs), the predicted Ha flux for the MS is significantly lower than what is found observationally 
(Weiner et al. 2002; Bland-Hawthorn & Maloney 2002), suggesting that an additional source of 
ionization must be invoked. Several possible mechanisms have been suggested, including an inter- 
action between the MS gas and the hot Galactic Halo gas, turbulent mixing, the existence of young 
massive stars embedded in the MS (Bland-Hawthorn & Putman 2001), shocks, self-interaction of 
the MS gas, and magnetic fields (Konz et al. 2001). 

The MS appears to be the result of interaction between the Galaxy and the MCs, but there 
is no consensus as to the exact form of this interaction (Putman 2000b). Theories have swung 
back and forth on the relative importance of tidal stripping (Murai & Fujimoto 1980; Gardiner 
& Noguchi 1996) and various kinds of gas dynamical interactions (Mathewson et al. 1987). The 
recent discovery by Putman et al. (1998) of a counter HI stream, leading the MCs, lends weight to 
the tidal origin theory. 

The work presented in this paper was motivated by the following two questions. 

(1) Is the MS dissipating into the Halo? 

The most fundamental issue about the origin and structure of the MS is to what extent interaction 
with the Galactic Halo determines or influences the MS gas. This problem becomes particularly 
important at the extreme northern end of the MS, because without pressure from an external 
medium the MS clouds should dissipate on a short time scale. However, if the external pressure 
is sufficient to confine the MS clouds, then the gaseous Galactic Halo has a much higher scale 
height than most theories predict. This paper addresses several issues in relation to this particular 
question. 

(2) Is a hierarchy of structures present in this almost primordial environment? 

In the Galaxy, a hierarchy of structure is present in the diffuse interstellar medium (ISM) , down to 
very small scales. This is mainly ascribed to interstellar (IS) turbulence. This process is seen both 
in the ionized medium, traced by pulsar scintillation and dispersion variations, and in the atomic 
gas, traced by the spatial power spectrum of HI emission and absorption fluctuations (Fiedler 
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et al. 1994; Green 1993; Lazarian & Pogosyan 2000; Dickey et al. 2001). Recently it has become 
possible to extend these studies to other galaxies, at least to our neighbors (Westpfahl et al. 1999; 
Stanimirovic et al. 1999; Elmegreen et al. 2001). The results suggest that the processes we see at 
work in the solar neighborhood are quite general. As the IS turbulence is driven, at least partially, 
by stellar processes such as winds and supernova remnants, it would be quite a surprise if the same 
turbulence spectrum was seen in an environment without any stars at all, like the MS. We have 
addressed this question briefly in Stanimirovic et al. (2002) and will concentrate on it further in a 
subsequent paper. 

This paper is organized as follows: In Section 2 we summarize previous HI observations of 
the MS, and especially its northern part. As parameters such as the MS age and its distance 
come up often throughout this paper, a brief description of various theories for the MS formation 
is also outlined here. Section 3 describes new HI observations of two regions at the tip of the 
MS conducted with the Arecibo telescope. The HI data and observational results obtained are 
presented in Section 4. Several clumps and their properties are shown in Section 5. In Section 6 
we discuss clump confinement issues and the constraints they place on the density of the Galactic 
Halo. Expected Halo densities, as well as several theoretical approaches for considering interactions 
between an MS cloud and the ambient medium, are discussed. An alternative explanation for the 
velocity field of one of the regions is also presented. We summarize our results in Section 7. 

2. Previous HI observations and theoretical models 

2.1. Lower resolution observations 

First HI detections of the MS were obtained by Dieter (1965) and Wannier &; Wrixon (1972), 
however Mathewson et al. (1974) were the first to discover a coherent large-scale structure which 
they named the Magellanic Stream. Observations with better resolution and sensitivity followed by 
Cohen & Davies (1975), Mirabel et al. (1979), Erkes et al. (1980), Mirabel (1981) and Mathewson & 
Ford (1984). These observations revealed the large-scale morphology of the MS, that is traditionally 
viewed as a long filament, containing a bead-like sequence of six discrete clouds: MS I, near the 
MCs at Dec -60°, to MS VI at the very northern tip at Dec ~ +15° (Mathewson et al. 1977), 
see the top panel in Fig. 1. Mirabel (1981) in particular focused on the tip of the MS, searching 
for evidence of its disintegration. Wayte (1989) used the Parkes 64-m dish to uncover complex and 
chaotic motions right at the MS tip. 

Fig. 1. — [fl.jpeg] The HI distribution of the MS from Mathewson &: Ford (1984) is shown in the 
upper panel. The lower panel shows clouds MS HI through MS VI of the MS as seen by the Parkes 
HIPASS survey (Putman et al. 2002). The two regions observed with the Arecibo telescope and 
presented in this paper are outlined with dashed lines. 
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The most recent HI observations of the MS, with spatial resolution of 15 arcmin and velocity 
resolution of 26 km s~^, were undertaken with the Parkes Multibeam system (Putman 2000b; 
Putman et al. 2002). These, as well as the recent observations obtained with the same instrument 
but with higher velocity resolution of 1 km s~^ by Briins et al. (2000), revolutionized the classical 
picture of the MS. Instead of six concentrations of gas, two striking, long and distinct filaments 
are seen that run together for most of the length of the MS but merge in several places, giving the 
appearance of a double-helix structure. North of Dec ~ 0°, the two filaments break into a network 
of clumps and smaller filaments (see the bottom panel in Fig. 1). A global velocity gradient is 
observed along the MS, from ~ +350 km s'^ (LSR) near the MC's to ~ -450 km s'^ (LSR) at 
the very northern tip. 



2.2. Theoretical models 

The two leading hypotheses for explaining the origin of the MS are based on tidal interaction 
versus ram pressure stripping, and are applied in several models. These two families of models 
predict different ages and distances for the MS. As the predicted parameters are important for the 
calculations in Sections 5 and 6 we carry out calculations for both sets of model predictions. 

The tidal hypothesis (Murai Sz Fujimoto 1980; Lin &; Lynden-Bell 1982; Gardiner et al. 
1994; Gardiner &: Noguchi 1996; Yoshizawa &: Noguchi 1999) — invokes the idea of gravitational 
stripping of the SMC gas by the Galaxy. According to N-body simulations, the MS was formed 1.5 
Gyr ago in a tidal encounter between the LMC and the SMC at their perigalactic passage. This 
encounter drew gas out of the SMC that later evolved into a leading bridge and a trailing tail. 
The MS follows the MCs on an almost polar orbit with a transverse velocity of 220 km s~^ (LSR). 
Simulations show that the MS gas occupies a range of distances, from 45 kpc in the vicinity of MS 
III, to 60 - 70 kpc at MS VI. 

One of the main drawbacks of the tidal model used to be the observation of only one tail of gas 
behind the MCs. Recent Parkes observations (Putman et al. 1998) have discovered the 'Leading 

Arm Feature' (LAF), a counter-stream that leads the direction of motion of the MCs, as expected in 
the tidal models. The absence of stars in the MS is another troubling issue for tidal models, as both 
stars and gas should be affected equally by tides. However, the most recent N-body simulations 
by (Yoshizawa & Noguchi 1999) show that a very compact initial configuration of the SMC stellar 
disk could result in only gas being disrupted, while the stars are left unaffected. 

The ram-pressure hypothesis (Moore k Davis 1994; Heller & Rohlfs 1994; Sofue 1994) 
— invokes the idea of ram-pressure stripping of gas from the MCs by an extended halo of diffuse 
gas around the Galaxy. Here, the MCs entered the extended ionized halo of the Galaxy some 500 
Myr ago at a galactocentric distance of 65 kpc. Material was then stripped from the SMC and 
the intercloud region and began to fall toward the Galaxy. In this model, the gas with the lowest 
column density, that lost the most orbital angular momentum, has fallen the furthest toward the 
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Galaxy. This gas is currently located at the tip of the MS, at a galactocentric radius of 25 kpc. 
The model requires that the diffuse Galactic Halo has densities of ~ 5 X 10""^ cm~^ at distance of 
65 kpc in order to strip gas, but not decelerate it too much. 

The tidal + drag model (Gardiner 1999; Moore & Davis, in preparation) The recent ob- 
servational confirmation of the LAF provides a new set of parameters for theoretical models. In 
order to better constrain the observational characteristics of the LAF, Gardiner (1999) introduced 
a weak non-gravitational term (a drag term) into the classical tidal model. The new simulations 
show interesting results: while the LAF is reproduced, simulated particles form a complete ring 
from the LAF, connecting the LAF to two trailing streams. These streams have different galacto- 
centric velocities starting from MS I to MS V, and then meet together again around MS VI (see 
Fig. 4 from Gardiner 1999). This simulation predicts that large-scale velocity splitting should 
be seen throughout most of the MS. The success of this model, based on both gravitational and 
non-gravitational forces, encourages further developments in the field that are eagerly anticipated 
(Moore, in preparation; Maddison, Kawata, &; Gibson 2002). 



3. Observations and data reduction 

The observations were obtained with the 305-m Arecibo telescope^ in June and July of 2000 
and June of 2001. The Gregorian feed was used with the narrow 21-cm (L-band) receiver. The 
illuminated part of the 305-m dish covers an area of about 210x240 m, giving a beam FWHM of 
approximately 3M x 3'. 4 arcmin (Heiles et al. 2001). Two observing bandwidths of 12.5 and 6.25 
MHz were used simultaneously on two correlator boards, each with 1024 channels, for each circular 
polarization. The resultant velocity resolution is 2.6 and 1.3 km s~^, respectively. 

Data were taken in the on-the-fly mapping mode, by driving in RA, at a rate of 4 arcsec per 
sec, and stepping in Dec by 2 arcmin between strips. The correlator was set to make one scan per 
strip by recording data continuously every 30 sec, which is effectively every 2 arcmin on the sky. 
As a result the data were slightly undersampled. At the beginning of each scan, a noise diode of 
known temperature was fired for 3 sec for calibration. 

Two regions of interest were selected from Mirabel et al. (1979): one centered at the northern 
tip of the MS, {l,b) = (86.7, -42.8), and the other in the MS V region, {l,b) = (94.4, -50.6) (Fig. 
1, these regions are outlined in the lower panel). Below we refer to these BifGclS clS MS VI and MS V, 
with the understanding that our maps do not cover the full extents of those clouds. Both regions 
were mosaiced with many small, overlapping maps, each of 10x10 arcmin, which take ~ 1 h to 
observe. The mosaic pieces were combined during the gridding process. 

We removed a linear baseline by fitting ranges of emission-free channels on both sides of the 



^The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated by Cornell 
University under a cooperative agreement with the National Science Foundation. 
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line. The gain correction, as a function of both azimuth and zenith angle, was then apphed. All 
spectra were subsequently convolved onto a sky grid using a Gaussian convolution function with 
FWHM of 2.6 arcmin, broadening the angular resolution to 4 arcmin. A gridding correction factor 
of 1.45 (Sramek & Schwab 1989) was applied to account for the new beam size (S. Stanimirovic 
2002, in preparation). 

The relationship between flux density and brightness temperature is given by ^(Jy beam~^) = 
O.lirb(K). The noise level is 0.04 K per 1.3 km wide channel, which is equivalent to a column 
density of ~ 10^^ cm~^. The peak brightness temperature is 1.5 K in the MS VI region, and 1.1 
K in MS V. For comparison, Mirabel et al. (1979) obtained a maximum column density of ^ 
4x10^^ cm~^ in the MS VI region using the Jodrell Bank telescope. In the same region we derive 
a maximum column density of 4.3x10^^ cm~^. 

4. Observational Results 

4.1. Morphology and velocity field of the MS VI region 

The MS VI data is shown in Fig. 2 where several velocity channels through the spectral-line 
data cube are presented. Fig. 3 shows position-velocity cuts. Fig. 4 shows maps of the HI column 

density and the intensity- weighted velocity field (panels a and b) . This region has several elongated 
clumps, inter-connected and forming a long curved filament. Often, sub-clumps are found nested 
inside larger clumps. We identify individual clumps and describe their properties in Section 5.1. 
Several interesting features appear in the data. 

A large loop-like feature, centered at RA 23"^ 07™ 26^ Dec 12°47' 30", with radius of 50 arcmin, 
is visible between i'lsr = —360 and wlsr = —385 km s"-*^. The loop is present through ~15 velocity 
channels but its size appears to stay the same from channel to channel. Although at first glance 
this loop is reminiscent of HI holes associated with expanding shells of gas seen in many galaxies, 
because no systematic change of the radius with velocity is found this cannot be interpreted as an 
expanding shell. Velocity profiles along the loop are broad but always single peaked. 

It is interesting to note that the most prominent clumps have a dominant orientation, from 
south-west to north-east, with position angle of ~ —40° (measured counter-clock wise from north). 
This corresponds to a position angle of ~ 90° in the Magellanic coordinate system, defined by 
Wakker (2001) and is parallel to the main axis of the MS (Putman 2000a). Similar orientation of 
clumps in MS VI was noted previously by Mirabel et al. (1979). 

The total HI mass in the mapped region is 3.7x 10^ Mq, if the distance of 60 kpc is assumed. 

Fig. 2. — [^.jpegjThe RA-Dec images of the MS VI region for different LSR velocities given in 
the top-left corner. The grey-scale intensity range is —0.05 to 0.3 K with a linear transfer function. 
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and it comprises about 28% of the HI mass in the whole MS VI region estimated by Putman 
(2000a). 

The velocity profiles in this region are single-peaked but quite broad with the dispersion being 

~ 15 km s~^. Several RA- velocity diagrams are shown in Fig. 3. At about Dec 13°05' an interesting 
velocity gradient of almost 40 km s~^ (LSR) is obvious. This gradient stays almost the same in 
the galactocentric velocity system, meaning that a spatial variation of the projection of the LSR 
velocity relative to the Galactic Center cannot explain the velocity field in this relatively small, 
40x40 arcmin, field. 

The velocity field shown in Fig. 4 is the first moment map of the HI cube. The two prominent 
features in this map are north-south ridge-lines at velocity —365 to —375 km s~^ in the east, and 
velocity —395 to —400 km in the west. The points labeled A and B in Fig. 4 indicate the 
extent of the east ridge which may continue beyond the edge of our map. Points labeled C, D and 
E indicate the extent of the west ridge. The steep velocity gradient evident in the middle panels 
of Fig. 3 is in the region bounded by points A, B, C and D in Fig. 4. This region alone resembles 
the velocity pattern of a rotating disk, however the velocity field of the larger area does not fit this 
interpretation. We discuss this further in Section 6.6. 

The velocity field in Fig. 4 is most likely caused by a combination of blending of smaller clumps 
along the line of sight and an interaction between the MS clouds and the ambient medium. With 
the Arecibo beam we are able to resolve several clumps blended together and forming a velocity 
gradient, see Fig. 3, especially around Dec 13° 08'. Some of these clumps show tails in velocity 
and sharp edges. This may be due to interaction with the surrounding medium. We discuss clump 
morphology in Section 6.7. Mirabel et al. (1979) saw a similar phenomenon in another MS VI 
region, while Wayte (1989) attributed broad velocity profiles to the breaking of the MS clouds into 
sub-clouds with different velocities. 

4.2. Morphology and velocity field of the MS V region 

Velocity channel maps presented in Fig. 5 show the general morphology of the MS V region. 
Several dominant large-scale features are visible. In the —320 km plane there is a round cloud 
centered around RA 23'' 39™, Dec 07° 35' 00". This region shows single-peaked velocity profiles. A 
long, straight filament stretches north-south for almost one degree, at RA 23'^ 38°' 30" from Dec 
8° to 9°, at velocity of ~ —327 km s~^. At higher velocities of —340 to —355 km s^^, the most 
prominent feature is a curved filament, centered at RA 23^^ 39°^ 30^, Dec 07° 40'. This consists of a 

Fig. 3. — [f3.jpeg] Position-velocity diagrams of the MS VI region taken at different Dec values, 
given in the top-left corner. The grey-scale intensity range is 0.05 to 0.7 K with a linear transfer 
function. 
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core superposed on an arc that extends for almost 1.3°. 

North of Dec = 07°20' all line profiles show double peaks. This is shown in the position-velocity 
diagrams, Fig. 6, where a single velocity component is present only in the upper left panel. All 
other panels show double profiles. The two velocity peaks are roughly centered at ~ —320 and 
~ —360 km s~^ and run almost in parallel. 

Double-peaked velocity structure has been seen previously in MS V by Mirabel et al. (1979) at 
very similar central velocities. Waytc (1989) also noted this velocity bifurcation in several regions 
in the MS. The most recent HIPASS data, with very coarse velocity resolution of about 26 km s^"^, 
show that spatially different filaments are common throughout the whole MS and most likely have 
different origins (Putman 2000a) . Similar line-splitting is found in several compact HVCs by Braun 
k, Burton (2000), who suggest that this structure is caused by organized outflows or inflows. 

What causes the dual velocity structure seen in MS V? The velocity separation between the 
two peaks of ^ 40 km s~^ is reminiscent of velocity features seen in the SMC. This was formerly 
interpreted as a signature of the SMC being ripped into two clouds which began during its last close 
encounter with the LMC some 200 Myr ago (Mathewson & Ford 1984). Recently, high spatial and 
velocity resolution observations of the SMC (Stavclcy-Smith ct al. 1997; Stanimirovic ct al. 1999) 
discovered more complex velocity profiles and numerous expanding shells of gas. This line-splitting 
in the SMC is now attributed primarily to the effect of expanding shells. Multiple structures in 
velocity are also seen in the Magellanic Bridge region (McGee & Newton 1986). No evidence for 
shells has been found in the MS. 

We have fitted overlapping Gaussians to the blended velocity components, and the results show 
that their velocity separation is quite constant over the area of MS V, as seen in Fig. 6. 

4-2.1. Gaussian decomposition 

To model the velocity profiles with one or a superposition of two Gaussian functions, the data 
cube was binned spatially in squares of 5 by 5 pixels. As an example. Fig. 7 shows one velocity 

Fig. 4. — [f4.jpeg] (a) An HI column density image of MS VI overlaid with contours. The grey-scale 
range is to 2.4x10^^ cm^^. Contour levels are from 7 to 30x10^^ cm^^, with a contour step of 
6x10^^ cm-2. The maximum column density is 3.9x 10^9 cm^^ at RA 23^' 07™ 52% Dec 11° 31' 
28". (b) The first moment map of the MS VI region overlaid with contours. Contours range from 
—400 to —360 km s~^, with a step of 10 km s~^. The grey-scale range is —391 to —364 km s~^ with 
a linear transfer function. Two prominent velocity ridges are shown with dashed lines: the east 
ridge with a mean velocity of —397 km s~^, indicated with points labeled with A and B; the west 
velocity ridge with a mean velocity of —370 km s~^, indicated with points labeled with C, D and 
E. An increase in noise toward map edges is noticeable. 
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Fig. 5. — [f5.jpeg]The RA-Dec images of the MS V region for different LSR velocities given in the 
top-left corner. The grey-scale intensity range is —0.01 to 0.9 K with a linear transfer function. 



Fig. 6. — [f6.jpeg]Position-velocity diagrams of the MS V region taken at different Dec values, 
given in the top-left corner. The grey-scale intensity range is to 1.3 K with a linear transfer 
function. 
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Fig. 7.— An example of Gaussian decomposition at RA 23^ 39°" 46", Dec 08°01' 00". Two different 
fitted components are shown as dashed and dot-dashed lines, while their sum is shown as a solid 
line. 
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profile, its decomposition into two Gaussian functions, and the residuals after subtracting the fit 
from the data. The mean rms vahic of the residuals in the whole region is close to the noise level, 
with its maximum below the 2-0" level. The two separate components arc centered around —350 
km s~^, the 'higher' velocity component ('H'), and —320 km s~^, the 'lower' velocity component 
('L'). The integrated intensity for each velocity component is shown in Fig. 8. The higher velocity 
component is the core-arc feature described above, plus a diffuse feature to the north-west. The 
lower velocity component is the round cloud and the north-south filament. The HI mass of the 
higher velocity component is 2.3x10^ M0. The HI mass of the lower velocity component is 3.2x10^ 
Mq (if the distance of 60 kpc is assumed). These together make up only 15% of the HI mass of the 
entire MS V region, 3.9x10^ Mq (Putman 2000a). 

A slight gradient in velocity of about 10 km is seen from the north to the south and is shown 
by the slight tilt of the ridge lines in Fig. 6. These trace roughly ~ —360 to ~ —350 km for the 
higher velocity component, and ~ —325 to ~ —315 km for the lower velocity component, from 
right to left on the lower middle panels of Fig. 6. The dispersion of the lower velocity component 
has everywhere a value of about 10 km s~^. The dispersion of the higher velocity component shows 
a significant variation with position. In the north, the velocity profiles are very narrow, with typical 
velocity dispersion of about 6 km s~^. The dispersion increases towards the south reaching about 
13 km along the core-arc feature. 

We have experimented with fitting more than two velocity components to the data, however 
the resultant fits did not improve the residuals significantly. Neither MS V nor MS VI regions 
show narrow line components indicative of cold HI (with temperature of 30 - 300 K). This is in 
agreement with theoretical predictions by Wolfire et al. (1995), that no cold cores should be found 
in clouds in the Galactic Halo at heights above 20 kpc. 

5. Clump analyses 

5.1. Clump identification 

In order to study the properties of the MS clumps and the indirect clues they can provide 
about the properties of the Galactic Halo, we have attempted to identify individual clumps. We 
used two different approaches. 

(i) To identify smaller clumps and sub-clumps within larger clumps, the three-dimensional 
clump-finding algorithm, 'clfind' coded in IDL (Williams et al. 1994) was applied to both data 
sets. A spatial binning by 4 pixels, and a velocity binning by 2 pixels was applied to reduce 
the computing time. 'Clfind' views data as a set of contour levels. It locates local maxima and 
defines clumps by tracing them by connecting pixels, at each contour level, that are within one 
resolution element of each other. The algorithm does not assume any particular geometry for the 
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Fig. 8. — An integrated intensity image of two Gaussian components fitted to MS V velocity 
profiles. The grey-scale intensity range is 1 to 4 xlO^^ cm~^ for the higher velocity component, 
and 1.2 to 4 xlO^^ cm~^ for the lower velocity component. 



Fig. 9. — [f9.jpeg]Three fully mapped and well isolated clumps in the MS V region. For each 
clump, the left-hand side shows an HI intensity image overlaid with contours of the total column 
density, while the right-hand side shows a velocity profile at the clump's center. 



Fig. 10. — flO.jpeg] Four fully mapped and well isolated clumps in the MS VI region, from the 
top to the bottom. For each clump, the left-hand side shows an HI intensity image overlaid with 
contours of the total column density, while the right-hand side shows a velocity profile at the 
clump's center. 
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clumps. The most important input parameter for 'clfind' is the contour increment. This was set to 
twice the noise level, as recommended in Williams et al. (1994). All selected clumps were further 
vigorously inspected by eye by overlaying the 'clfind' output on the original data using the KARMA 
visualization package (Gooch 1995). Only a few (seven) clumps that are isolated and fully inside 
the edges of our map were considered for further analysis. Those clumps, three in MS V and four 
in MS VI, are shown in Figs. 9 and 10. For each clump the central position, central velocity, clump 
size (i?c) measured from the integrated intensity map, FWHM velocity line width {Av), and peak 
brightness temperature were estimated and are listed in Table 1. 

(ii) The clumps found by 'clfind' may or may not be distinct physical structures. On larger 
scales, all of MS VI and the 'L' and 'H' components of MS V could be considered as clumps 
themselves. The corresponding quantities for these three largest clumps are also listed in Table 1. 

5.2. Observed and derived properties 

Fig. 11 shows the angular size and line width for all the clumps. A correlation between line 
width and cloud size has been found for interstellar clouds in a number of surveys (references in 
Vazquez-Semadeni et al. 1997). This is one of the famous Larson scaling laws: Av oc Re- For the 
case of self-gravitating molecular clouds (3 = 0.5 (Myers Sz Goodman 1988), a slightly flatter slope 
of 0.4 was found by Falgarone et al. (1992) in molecular clouds that are very little contaminated 
by star formation. We overlay the Larson relation in Fig. 11 for the assumed distance of 60 kpc. 
While all sub-clumps and one of the largest clumps (the whole of MS VI) appear to follow the 
correlation, two other large clumps have Av lower by about 10 km s^^ than the predicted values. 

We now derive several important clump properties. These properties depend on the assumed 
distance of the MS. Figs. 12 through 16 have double axes corresponding to the extreme distances 
of 60 and 20 kpc (see Section 2.2). For all these figures, the left-hand and bottom axes correspond 
to the far distance, and the right-hand and top axes correspond to the near distance. 

The total HI mass is derived by integrating the brightness temperature over the entire line 
width. The HI mass is plotted as a function of cloud size in Fig. 12. For a given size our noise 
level sets a minimum detectable cloud mass. Using the narrowest measured clump line widths of 

~ 10 km s~^ we get the threshold indicated by the dashed line on Fig. 12. Fig. 12 suggests that 
there may be many more structures like our clumps in MS V and MS VI which could be found by 
a more sensitive survey. 

Assuming that clouds arc spherically symmetric, using their HI mass and radius we can es- 
timate an average HI volume density, ric- We find the values of shown in Fig. 13, where the 
3-0" detection threshold is again shown by the dashed line. A typical value for the small clumps 
is 0.05 cm~^ and for the large clumps 0.02 cm~^, for the 60 kpc distance assumption. For the 20 
kpc distance the volume density increases by a factor of three, as shown by the right-hand scale on 
Fig. 13. 



-13- 





40 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

^MS VI sub-clumps 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


T 


35 


V sub-clumps 




^ M S/¥l - 


C/) 








"WHM, Av (km 


30 
25 




- 

^MS VH \ 


lump ' 


20 


: O 




V MS VL : 


o 


15 


/ill 1 1 1 1 




II 






10 20 


30 


40 50 






Clump Angu 


or Radius, 


(arcmin) 
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6. Discussion 



6.1. Cloud confinement 

We now turn to cloud stability issues and the survival of the MS so far away from the MCs. 
There are several questions to ask about the HI structures which we find. 

(1) Are the MS clumps gravitationally confined? 

In order to be gravitationally confined, a clump's escape velocity must be larger than /S.v/2. This re- 
quires that the clump's total mass is Mgrav > RAv^/8G. Fig. 14 compares the measured HI masses 
of our clumps with this hypothetical gravitational mass. For all our points 30 < Mgj-av/^m < 250, 
if the distance of 60 kpc is assumed, and 120 < Mgrav/-^Hi < 750, for the 20 kpc distance assump- 
tion. Thus gravitational confinement would require unreasonable amounts of dark matter. 

(2) Are the MS clumps unbound and freely expanding? 

In the tidal model the MS age is 1.5 Gyr, while in the ram pressure model the age is 500 Myr. If 
the clouds are freely expanding their expansion time is given by: iexp = 2R/Av. These numbers 
are much less than the MS age. Fig. 15 plots the expansion age as a function of clump radius. 
The expansion age as a fraction of the MS age is indicated by the horizontal lines which apply to 
either model. This huge discrepancy between the expansion time and the age is the fundamental 
problem of the gas dynamics of the MS, as pointed out first by Mirabel et al. (1979). Barring the 
unreasonable dark matter halos required for gravitational confinement the simplest explanation is 
that the MS clouds are confined by a hot gaseous Galactic Halo. 



6.2. Density of the Milky Way Halo 

If the MS is confined by the Galactic Halo, then the external pressure, = ^^h^li, must equal 

the clouds' internal pressure Pc = kucTc'. Pc < Ph- The relevant value of is not the microscopic 
kinetic temperature but the total random motion of the atoms which we describe as a temperature 
defined by the line-width: 

k f Av V 

Tc = — — = 21.9 K X ^ . 1 

81n2mH Vkm s'^ J ^ ^ 

Assuming a halo temperature of 10^ K and using Av for the clumps in Table 2 and the densities 
Uc illustrated on Fig. 13, we get Halo densities shown in Fig. 16 as a function of distance from 
the Galactic Plane, z. Typical values for the near distance {z 15 kpc) are = 10~^ cm~^ and 
for the far distance {z ~ 45 kpc) are rih = 3 x 10"^ cm~^. Note that these estimates are upper 
limits on the Halo density as magnetic, turbulent or ram pressure contributions were not included. 



-16- 



0.08 



o 

Q. 

O 
CD 



0.06 - 



0.04 



o 



0.02 



1.8 



log(/?.) (pc) [ d = 20 kpc ] 
2.0 2.2 



2.4 



1 1 1 

• 






0MS VI - 

Qus V " 


0.25 










0.20 


• 










\ 


o ° 






0.15 
0.10 


1 , , 


'iniif 

,1,1, 


MS VI jx 
MS VH ^ I 

MS VL y 

1 ~ 7 ~ ^ _ , T 


0.05 



2.2 



2.4 2.6 2.8 

log(/?,) (pc) [ d = 60 kpc ] 



o 

O 
(N 

II 



O 



3.0 
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6.3. Expected Halo density 

Current theories of galaxy formation predict the existence of a hot Galactic Halo containing 
the leftover gas which was unable to cool since the formation of the Galaxy. However, the density 
structure of the Halo is still very uncertain, especially at large distances from the Galactic Plane. 
Here we compare our results with several previous estimates, both observational and theoretical. 

6.3.1. Observational constraints 

The only direct observations of the hot Halo gas come from measurements of the soft X-ray 
background (Snowden et al. 1997). These X-ray maps reveal hot gas with a temperature ~ 10^'^ 
K, a scale-height of 1.9 kpc and an electron density of ~ 3.5x10^'^ cm^'^ Indirect observational 
constraints on the Halo density are more common. These are based on dispersion measure mea- 
surements of pulsars in globular clusters, or on the observational properties of HVCs, including the 
MS clouds. 

Mirabel et al. (1979) estimate that nh = 7 x 10~^ - 7 x 10~^ cm~^ for the distance range 10 - 
100 kpc from HI observations of the MS and under similar confinement considerations as discussed 
in this paper. Similar vahics, nji > 10~^ cm""', were suggested from Hq observations along the 
MS by Weiner & Williams (1996). Quilis &: Moore (2001) applied the terminal velocity model by 
Benjamin &: Danly (1997) for HVCs moving at ~ 300 km s~^, at distances 10 - 100 kpc and with 
HI column density of lO^'' cm~^. This results in the requirement that nh(10 kpc) < 3 x 10^^ cm~^ 
and nh(100 kpc) < 3 x 10~^ cm~^, under the assumption that these are pure gas clouds. Our 
results agree well with the estimates from these three studies. 

Another constraint for the Halo density at the distance of 130 kpc of ~ 10~^ cm~^ comes 
from HI observations of HVC 125-1-41—207 by Briins et al. (2001). A significantly lower value for 
the upper limit on the Halo density, nh(50 kpc) < 10~^ cm~^, was imposed by Murali (2000) by 
allowing the whole MS VI cloud to survive at the tip of the MS for about 500 Myr. We look more 
closely into this approach in Section 6.4. 

6.3.2. Theoretical predictions 

Benjamin & Danly (1997) considers a Galactic Halo consisting of three components: the warm 
ionized layer of HII parameterized by Reynolds (1993); the mean HI density distribution given by 
Dickey & Lockman (1990); and the hot Halo component with Th = 10^ K prescribed by Wolfire et al. 
(1995). While the first two components were parameterized observationally, the third component 
was constructed theoretically assuming an idealized isothermal halo in hydrostatic equilibrium, with 
the mid-plane density matching X-ray emission data by Garmire et al. (1992). The combination 
of all three components results in the analytic representation for as a function of z given by 
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Benjamin & Danly (1997): 



nh(z) = 1.1 X 10-'^ [1+ (4p,/19.6)] 



-1.35 



-3 



(2) 



cm 



At z = 15 kpc the expected density is nh = 4 x 10 ^ cm and at z = 45 kpc nh = 2 x 10 ^ cm ^. 
These values are 25 - 150 times lower than those found in Section 6.2. 

The ram-pressure model for the MS origin by Moore & Davis (1994) requires slightly higher 
values for the Halo density in order to match the kinematics of MS clouds: 



Hence, at a distance of 20 kpc this yields a gaseous halo with nh = 3 x 10~^ cm~^, while at a 
distance of 60 kpc nh = 6 x 10~^ cm~^. Our estimates are much closer to the values predicted by 
this model but still several times higher. 



There are (at least) three different approaches present in the literature concerning the dynamics 
of clouds moving through an ambient medium. As understanding of the nature of this interaction 
can in turn provide clues about the properties of the Halo, here we discuss these approaches briefly. 

(1) The mass loading approach 
As the MS moves through the hot ambient Halo gas its clouds are subject to a drag force. Murali 
(2000) showed that this motion is actually dominated not by drag but by strong heating from 
accretion. If the Halo density is too high, heating due to accretion can cause cloud evaporation. 
This process is usually referred to as mass-loaded flow and it is common in comet interactions with 
the solar wind (Flammer 1991). The Halo gas ionizes the cloud's leading edges and picks up the 
'new material' in a very asymmetric way that is constrained by the magnetic field: magnetic field 
lines slip past the cloud therefore allowing ablation of material only from the poles. Murali (2000) 
accounted for these effects and estimated that for the survival of the whole MS IV region for 500 
Myr, at a distance of 50 kpc, nh < 10~^ cm~^ is required. The key factor in this approach is the 
cooling resulting from a cloud's mass-loss. 

Following Murali's calculations (their equation 1) for the MS V and MS VI clouds studied here, 
drag and accretion set a very low upper limit for the Halo gas density: nh ~ 10~^ cm~^ at z = 45 
kpc, and nh ~ 10~^ cm~^ at z = 15 kpc, if mass-loss is not included. These values are significantly 
lower than theoretical predictions. Unless a significant cooling due to mass-loss is introduced, it is 
very hard for the hot Halo gas to exist in this model. 

If the mass-loss rate, M= 2.5 x 10~^ M© yr~^ (Murali 2000) is included, the requirements are 
comparable to our estimates: nh ~ 3 x 10~^ cm~^ at z = 45 kpc and nh ~ 3 x 10~^ cm~^ at z = 15 



nh = 0.002/ [1 + {R/12 kpc)f cm'^. 



(3) 



6.4. 



Nature of Galactic infall 
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kpc. A slightly lower rate of M= 10~^ Mq yr~^ results in ten times lower Halo densities. However, 
the mass-loss rates of M= 10~'^ or fewxlO^^ Mq yr^^ throughout the whole cloud life-time are too 
high in comparison with the cloud HI masses discussed in Section 5.2. Such mass-loss rates would 
imply that cloud HI masses were 10 - 1000 times higher at the time of the MS formation. Smaller 
rates that are in agreement with the current HI masses are not sufficient to allow the existence of 
the Halo gas. Detailed simulations of the mass-loaded mechanism may help to constrain parameters 
such as mass-loss rate. 

(2) Ballistic approach 

Gregori et al. (1999) showed results of a 3-D study of a ballistic interaction of a moderately super- 
sonic dense cloud with a warm magnetized medium. Since clouds are supersonic, their motion leads 
to the formation of a forward bow shock and a reverse crushing shock propagating through the 
cloud. As the cloud moves, its surface is subject to several instability mechanisms, the most dis- 
ruptive of these being the Rayleigh- Taylor instability (R-T) that develops at the interface between 
two flTiids when the lighter fluid accelerates the heavier one. These instabilities slowly disrupt the 
entire cloud and can dramatically change its morphology. In contrast to the mass-loaded flows, 
magnetic field lines here stay trapped, causing the development of strong magnetic pressure at the 
leading edge of the cloud. 

The time-scale over which this dramatic change occurs is several times the so called 'crushing 
time', TcR = 2RcX^/'^ /Mcs, with M being the Mach number, x being the ratio of cloud to ambient 
density, and Cg being the sound speed. During a time shorter than tcr, the magnetic pressure is 
comparable to the ram pressure, while afterwards ram pressure dominates. For the MS clouds, 
M ~ 2, and assuming x = 100 from Section 6.2, tcr ~ 10 - 30 Myr, depending on whether the 
near or far distance to the clouds is assumed. As a few x tcr is a significant fraction of the MS 
age this suggests that during a large fraction of a cloud's life-time magnetic pressure may play a 
significant role. Inclusion of magnetic pressure in the calculations in Section 6.2 may somewhat 
reduce the required values for the Halo density. 

(3) Hydrodynamical approach 

Full hydrodynamical, 3-D simulations of HVCs moving through a diffuse hot gaseous component 
were performed recently by Quilis &; Moore (2001). They investigated both pure gas clouds, being 
in pressure equilibrium with the external hot medium, and dark matter dominated clouds with 
an additional potential field that maintains their dynamical equilibrium. Particular cometary 
morphology, seen in the case of many HVCs as well as in many MS clouds, was reproduced for both 
cloud types but under the condition that the diffuse medium has density nh > 10~^ cm~^. 
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6.5. Terminal velocity and neutral gas fraction 

At the tip of the MS, clouds have had plenty of time to decelerate by a drag force (Murali 
2000). The present velocity of the MS clouds is moderately supersonic (220 km vs. sound speed 
Cs of ~ 100 km s~^) suggesting that clouds are most likely now at the terminal velocity (Benjamin 
and Danly 1997). This means that clouds have stopped decelerating and that their boundary layers 
are subject to Kclvin-Hclmholtz instabilities, while Rayleigh- Taylor instabilities are not important 
at this stage. Benjamin &: Danly (1997) show that for a population of clouds with known distances, 
terminal velocities and column densities, it is possible to constrain the mean density of the gaseous 
halo and cloud neutral fraction fc = iVHi/(A^Hi + -^Hli): 

2Nnig{z) 

C-ofcnh = 2"^' W 

here Cd is the drag coefficient, g{z) is the gravitational acceleration, and is the cloud's terminal 
velocity. The mean column density for the MS V and MS VI clouds is (logA^Hi) ~ 19 cm~^. 
Assuming that vt = 220 km s~^, the right-hand side of the equation does not depend on distance 
and has a mean value of {C-ofcPh) ~ 6 x 10~^ cm~^. If we assume Cd = 1 and apply Halo densities 
from Section 6.2, an estimate of the cloud's neutral fraction can be obtained: at z = 15 kpc fc = 0.6, 
while at 2; = 45 kpc /c > 1. 



6.6. An alternative, gravitational confinement 

An alternative explanation for cloud confinement is that they are gravitationally bound. This 
paradigm has been proposed for several compact HVCs (CHVCs) which are in some respects similar 
to MS V and MS VI clouds. An additional hint for this is provided by the peculiar velocity field 
seen in the MS VI region (Figs. 3 and 4), and discussed in Section 4.1, that may be interpreted 
as being due to differential rotation. In the case of three CHVCs, Braun &; Burton (2000) found 
evidence for such orbital motion. In order to compare the MS VI cloud with these CHVCs we 
have attempted to model the velocity field north of Dec 12°30' using the tilted ring algorithm. The 
derived rotation curve has a quick and almost linear rise to ~ 15 km s~^, and then flattens out 
with large uncertainties. 

This curve was used to obtain an estimate of the dynamical mass: Mdyn = 1.5 x 10'' M© or 
Mdyn = 5 X 10^ Mq, depending whether the far or near distance is assumed. This further implies 
dark-to-visible mass ratios of 50 to 180, respectively, after including a contribution by helium of 
40% by mass to the observed HI mass. Therefore, if interpreted as a rotating HI cloud, this would 
require a very dark matter halo. 

No model for the formation of the MS has ever predicted the existence of dark-matter domi- 
nated clouds in this gas filament. Even the parent galaxies, the MCs, are not rich in dark matter. 
It is therefore hard to believe that rotational support is realistic. This was the conclusion of Sec- 
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tion 6.1. However the derived rotation curve is very similar to the rotation curves of CHVCs derived 
by Braun & Burton (2000). Dynamical masses of these CHVCs are about 10 times higher than the 
dynamical mass of the MS VI region, while the dark-to- visible mass ratio of the MS VI region is 
a few times higher than that for CHVCs. As the velocity field of MS VI is most likely the result 
of clump blending along the line-of-sight, it is interesting to note how well this effect can mimic 
differential rotation. 

6.7. On the clump morphology 

Several clumps in the MS data cube show gradients in the brightness and velocity distribution, 
exhibiting a cometary or head-tail morphology. The best example is the south-eastern side of the 
MS V cloud at —340 to —355 km s~^, see Fig. 17. This bow-shock-like feature may be a strong 
compression front, centered at RA 23*^ 40™, Dec 07°34' with two tails swept back on either side. 
This is strongly suggestive of a supersonic interaction between the HI cloud and an external, low 
density, ionized medium. 

Similar cometary morphology is seen in many HVCs in surveys of Meyerdierks (1991) and 
Briins et al. (2000). Putman (2000a) found many such clouds in the MS, which she speculated 
may be in the process of evaporation. The 3-D hydrodynamical simulations of an HVC moving 
through a diffuse hot medium by Quilis & Moore (2001) show a bow shock in front of the cloud, gas 
compression at the front edge of the cloud, and visible or invisible tails behind the cloud, depending 
on the density of the ambient medium. The double-tail morphology seen here is more similar to the 
simulations of the interaction between a disk-like structure with the surrounding medium by Quilis 
et al. (2000). These simulations were the first successful attempt to incorporate complex turbulent 
and viscous stripping at the interface of the cold and hot gaseous components. They predict a 
compression front and two tails folding on both sides. The similarity of the bow-shock-like feature 
with the simulations by Quilis et al. (2000) demonstrates the importance of turbulent and viscous 
mixing at the boundary layers when considering cloud - ambient medium interactions and should 
be fully explored in the future. 

7. Conclusions 

In this paper, we have presented new HI observations of two regions, MS V and MS VI, at 
the northern tip of the MS using the Arecibo telescope. The new data sets show the complex 
morphology of the MS with numerous inter-connected clumps. We have found double velocity 
structures in the MS V region. This region also contains an interesting bow-shock-like feature that 
is strongly suggestive of an interaction between the MS and an external medium. The MS VI data 
set shows a large velocity gradient, caused most likely by blending of small clumps along the line 
of sight. This region resembles the velocity field of a rotating disk, but we find that rotation is not 
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Fig. 17. — An HI intensity image of the bow-shock-like feature in the MS V region at the LSR 
velocity of —351 km s~^. 
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a likely explanation. 

Several MS clumps have been isolated to investigate their confinement mechanism. We show 
that unreasonably large amounts of dark matter are required in order for clumps to be gravitation- 
ally confined. Clumps do not seem to be in free expansion cither. The easiest way to explain the 
clump properties is with external pressure confinement by the hot Galactic Halo. In this scenario 
we place an upper limit on the Halo density: nh = 10~^ cm""^ at z « 15 kpc, and nh = 3 x 10~^ 
cm~^ at 2 45 kpc. These values agree well with several previous indirect observations of the 
Halo, but are significantly higher than the values predicted theoretically for an isothermal hot Halo 
by Wolfire et al. (1995). Our results are closer to the expected Halo densities of the Moore k, Davis 
model. 

Although cloud evaporation is most likely an important process at the tip of the MS, the 

mass-loaded approach for cloud interaction with the ambient medium requires very high mass-loss 
rates to enable the existence of the hot Halo gas. The ballistic consideration of clouds in the Halo 
suggests that during a large fraction of the cloud lifetime magnetic pressure may have a significant 
role in cloud evolution. Inclusion of magnetic pressure would reduce somewhat the required values 
for the Halo density. 

Hydrodynamical modeling of cloud interaction with the ambient medium also requires high 
Halo densities to reproduce cometary features commonly seen in HI observations of HVCs. If we 
assume that MS clouds are at their terminal velocity, the cloud's neutral fraction can be estimated 
at z = 15 kpc /c = 0.6, while at z = 45 kpc /c > 1- 
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Table 1: A summary of clump properties. 



Clump ID 


RA (J2000) 


Dec (J2000) 


LSR Velocity 


Radius 


FWHM 


^pecik 




(lili:mm:ss) 


(dd:mm:ss) 


(km s"-*-) 


(arcmin) 


(km s~^) 


(K) 


Small clumps 


MS V: 














1 


23:40:23 


07:58:45 


-354.5 


13 


21.1 


0.8 


2 


23:38:55 


07:38:45 


-318.3 


18 


18.6 


0.7 


3 


23:38:07 


07:36:45 


-344.2 


15 


24.2 


0.7 


MS VI: 














1 


23:07:07 


11:55:45 


-371.1 


13 


20.7 


1.0 


2 


23:07:55 


13:03:45 


-363.3 


11 


16.7 


0.8 


3 


23:07:23 


11:37:45 


-373.7 


12 


16.3 


0.8 


4 


23:06:19 


12:23:45 


-391.8 


13 


16.4 


0.6 


Large clumps 


MS VH 


23:40 


08:10 


320 


64x28 


23.5 


0.9 


MS VL 


23:39 


08:19 


~ -350 


54x30 


20.6 


1.2 


MS VI 


23:07 


12:22 


~ -370 


65x31 


35.5 


1.7 
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